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Abstract 

In this article we examine a variety of quantitative models for describing 
archaeological networks, with particular emphasis on the maritime networks 
of the Aegean Middle Bronze Age. In particular, we discriminate between 
those gravitational networks that are most likely (maximum entropy) and 
most efficient (best cost/benefit outcomes). 

1 Introduction 

In archaeology, the primary sources of information are usually finds from specific 
sites. In order to understand the social, cultural and political context of these finds 
it is essential that we understand how sites relate. However, it often requires so much 
effort to obtain the physical information from a single site that there is a danger 
that we become 'site-centric'. Despite our best attempts at deducing relationships 
from the artefacts found at them, there is often little direct information about how 
sites interact. Quantitative modelling provides one response to this challenge. Good 
quantitative modelling can be insensitive to poor data, make assumptions and biases 
clear and debatable, and can allow us to provide possible answers to questions that 
could not be asked any other way. In the best case, such answers can be checked 
later from the archaeological record. 

One of the major problems that we attempt to address is to determine how the 
relationships between (archaeological) sites are conditioned by geographical space; 
to what extent does the exchange between them transcend their geographical con- 
straints? As the basis for our modelling we have found the language of complex 
networks particularly useful as a network is both a set of vertices — the sites — and 
a set of edges, or links — the relationships between them. Socio-spatial networks 
do not arise and develop arbitrarily. Their maintenance incurs explicit and implicit 
'costs', both in sustaining sites and in supporting links [5]. The common agency in 
the networks that we shall discuss below is the need to regulate or distribute these 
costs. We now encounter a conundrum. Despite these general features, the mod- 
elling of social networks requires very specific tailoring for individual case studies. 
This is contingent on the level of social organisation, the nature of the exchange, 
relevant distance scales, travel technology, etc. of the society being studied. Thus, 
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for exemplary purposes, we have found it useful to work with one particular data 
set, namely the Middle Bronze Age (MBA) to Late Bronze Age I (LB I) S.Aegean. 
Nonetheless, we shall try to be as general as possible within this framework. 

2 Networks 

As we have said, networks primarily consist of vertices (or nodes), representing 
agents or populations or resources, and edges (or links), which represent the ex- 
change between them, from the physical trade of goods to the transmission of ideas 
or culture. The MBA Aegean has several distinct features that lends itself to a 
network analysis. There is a relatively clear temporal delineation both at the begin- 
ning of the MBA and at the end of LB I. From the start of the MBA (c.2000 BC), 
society on Crete, often dubbed 'Minoan', sees a series of innovations, such as mon- 
umental architecture at palatial centres, new craft technologies and more sustained 
long-distance exchange. This last trait is most likely enabled by the innovation of 
the sail that seems to have occurred around this time, supplementing what was 
previously a paddle-based maritime technology. The end of our period of study, 
LB I, c.1450 BC, sees the destruction of the Minoan palaces and a swing in power 
and influence to the Greek mainland; this marks a useful delineation for analytical 
purposes. The physical boundaries of our focus for study are also relatively well 
defined, though long range trade beyond the Aegean was essential e.g. in providing 
tin, and questions regarding the relationship to Egypt remain important. Neverthe- 
less, as a first approximation, limiting our analysis to the Aegean region provides a 
good example of where modelling may be of particular use. There is a large amount 
of material from this era, with Akrotiri and Minoan palaces such as Knossos being 
rich sources. In this paper we will describe our network model, ariadne, and show 
how it is related to other approaches to archaeological modelling, using these MBA 
maritime exchange networks as a testing ground. In so doing we will go beyond the 
qualitative descriptions given in our previous work [T21 [L7] . 

2.1 Defining Our Vertices 

One of the main ways in which the nature of society impacts on the construction 
of networks is in determining what constitutes vertices which, in principle, could 
vary in scale from individuals/households to communities to islands. Ironically this 
also reflects the fact that most network analysis is also site-centric (although graph 
theory does provide some tools to get around this [T3l HI])- in general, island 
archipelagos are ideal for networks as the geography provides a natural definition 
for our vertices. Our working hypothesis is that the detailed behaviour of individuals 
can be subsumed in the behaviour of the 'island' community. In this regard, we have 
taken each vertex to correspond to an island, or a large part of an island or, in the 
case of N. Crete and other coastal areas, to an isolated centre of population. The 
latter sites effectively behave as islands because of the difficulty of land travel. This 
hypothesis is contingent on a maritime technology permitting long-distance travel 
that enables the island (or 'island-like' coastal site) to be thought of as the basic 
unit for 'trade' and has been discussed in detail elsewhere [T2], [T7] . However, once we 
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move away from individuals to communities, it is apparent that each vertex has two 
very different attributes; its resource base and the exploitation of these resources. 
The former, the carrying capacity of the vertex, is part of our input. The latter 
is represented by the 'size' of the community that inhabits it. That is part of our 
output. 

In many models there is no discrimination between carrying capacity and popu- 
lation (i.e. one is taken proportional to the other) and, often, the carrying capacities 
of the sites are set equal. This may well be appropriate for some cases such as pre- 
urban civilisations where land was occupied in small villages of roughly equal size. 
The work of Broodbank [4J on the Early Bronze Age in the Cyclades is a good exam- 
ple of this approach. There the potential cultivable land of each island was assessed 
and the density of sites on habitable land was chosen to be uniform. However we are 
interested in an era when significant differences emerge between the size of different 
sites. These differences reflect more than just local differences in resources and we 
are looking to understand how interactions with other sites may have supported 
very large sites. Thus it is important that site sizes are variable in our model, both 
in terms of the fixed inputs given and in terms of the outputs, since we expect a 
complex interdependence of interactions and site sizes. We have chosen a set of 
major known islands or coastal sites as representing the most important locations 
for vertices. See Figure 1, with details in Tabled] The choice is made on the basis 
of the archaeological record. It is sufficient for our purposes in Table 1 to classify 
the carrying capacities of the sites as 'small' (S), 'medium' (M) or 'large' (L). What 
this means quantitatively will be explained later. In so doing we are assuming that 
other locations are peripheral to the dynamics of the whole system or at least their 
effect is well represented by the dominant site in their region. 




Figure 1: Thirty nine important Aegean sites of the MBA (Middle Bronze Age). 
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1. 


Knossos 


(L) 


14. 


Kea 


(M) 


27. 


Mycenae 


(L) 


2. 


Malia 


(L) 


15. 


Karpathos 


(S) 


28. 


Ayios Stephanos 


(L) 


3. 


Phaistos 


(L) 


16. 


Rhodes 


(L) 


29. 


Lavrion 


(M) 


4. 


Kommos 


(M) 


17. 


Kos 


(M) 


30. 


Kasos 


(S) 


5. 


Ayia Triadha (L) 


18. 


Miletus 


(L) 


31. 


Kalymnos 


(S) 


6. 


Palaikastro 


(L) 


19. 


Iasos 


(M) 


32. 


Myndus 


(M) 


7. 


Zakros 


(M) 


20. 


Samos 


(M) 


33. 


Cesme 


(M) 


8. 


Gournia 


(L) 


21. 


Petras 


(L) 


34. 


All 1 

Akbuk 


(M) 


9. 


Chania 


(L) 


22. 


Rethymnon (L) 


35. 


Menelaion 


(S) 


10. 


Thera 


(M) 


23. 


Paroikia 


(M) 


36. 


Argos 


(M) 


11. 


Phylakopi 


(M) 


24. 


Amorgos 


(S) 


37. 


Lerna 


(M) 


12. 


Kastri 


(M) 


25. 


los 


(S) 


38. 


Asine 


(S) 


13. 


Naxos 


(L) 


26. 


Aegina 


(M) 


39. 


Eleusis 


(M) 



Table 1: Thirty- nine sites of the MBA Aegean. Note Knossos as site 1 and Akrotiri, 
subsequently destroyed in the eruption of Thera, as site 10. The S,M,L indicates a rough 
assignment of size — small, medium or large. 



As for notation, we denote vertices using lowercase mid Latin indices i,j, .... The 
fixed carrying capacities of sites will be denoted Si (habitable land available). In 
some models the total population of a site may differ from this as we will indicate. 

2.2 Defining our Edges 

The edges represent the interactions and can be seen within the framework of dif- 
ferent types of space e.g. artefact space [2S], [2H1 E2] • However we take the view that 
archaeology is rooted in geographical space and it is an unavoidable constraint on 
interactions. 

This means that one of the most important inputs to our model will be a table 
of distances between our sites. These may be simple Euclidean (as-the-crow-flies) 
distances. A more sophisticated approach will use estimates of typical journey times. 
This will be sensitive to the technology available and could require information on 
currents, typical winds, slopes, and even security. For instance in our examples we 
have analysed the layout of islands in the Aegean by hand to take promontories 
into account, but not the effects of winds and currents. Likewise for land travel we 
picked routes by hand that reflected current three-dimensional geography, arguing 
that the accuracy of a full least cost path analysis is unnecessary given the level 
of approximation. Finally in the examples shown here, we choose to penalise land 
travel over sea travel by a 'friction' coefficient of 3.0. Results are largely insensitive 
to reasonable values of the coefficient (larger than unity), primarily affecting the 
width of Crete. We will not look at this issue in any depth here, rather we will take 
it that these effective distances between sites are given. The issue we address in 
detail is how to model the actual interactions, given a framework of sites and their 
separations. 

For notation, if the attribute of an edge from i to j is A%j then the matrix A of 
these elements is an adjacency matrix for the network. When Aij defines the flow 
from site i to site j it will be denoted by Fij. Note that unless otherwise stated, 
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models exclude self-loops excluded, that is Fa = 0. The effective distance from site 
i to site j will be c?y. For all examples here, this is symmetric dij = dji, although as 
discussed above this is not always appropriate. 

3 Geographical Models 

Our emphasis, and that of the other main models that we shall discuss, is based on 
network analysis, as we have already indicated. However, before going into details we 
stress that there are alternative approaches that work more with 'zones of influence' 
than exchange directly. Although vertices are preserved, networks themselves are 
abandoned, and we begin by discussing such models briefly. 

3.1 Geographical Models without Networks 

There are several examples of models for archaeology, based on significant sites 
such as we have defined, yet without any explicit edges defined and so no network. 
If we treat all sites as equal we can construct partitions (or partial partitions) of 
geographical space into something like 'zones of influence'. This is a very old idea, 
reincarnated in its simplest form as simple Voronoi, or Thiessen, polygons, each of 
which contain all the points in space for which the single site at the centre of that 
polygon is the closest site. They can be used to indicate zones of control of each 
site, e.g. for Etruscan Cities using Euclidean distances [23J. The XTent model of 
Renfrew and Level [21] is a generalisation for cases with different site sizes. Here 
site % of size Si is deemed to 'dominate' site j of size Sj if tan(#) > dij/(Si — Sj) 
where d^ is the distance between the two sites and tan(0) is a parameter of the 
model. In the XTent model one can work with just the sites (rather than all points 
in space needed for Voronoi diagrams) and we can represent the resulting hierarchy 
of sites as a directed network. An example of how this may be used with modern GIS 
techniques to gain a good estimate of actual walking time, rather than using simple 
Euclidean distance, has been given for Neopalatial Crete by Bevan J3J. However, 
such a network representation adds little and the model is not usually visualised in 
this way, e.g. see [2~i] 13]. 

3.2 Simple Geographical Networks 

There are two simple ways to construct a network which captures non-trivial infor- 
mation about the global interactions between sites, both based on a thresholding of 
the distance matrix. 

The first is a 'maximum distance network' (MDN) in which an edge from site 
i to site j is connected if the distance is less than some model parameter D. In 
the language of links, we might construe this as imposing a link likelihood Aij = 
Q(D-dij) for the journey, where Q(x) is the step function, taking value 1 for dij < D 
and zero otherwise. MDN is used in many fields to establish distance scales, e.g. as 
a model of ad-hoc wireless networks formed between mobile devices [30]. A great 
deal of work exists on these types of models (particularly 'random geometric graphs' 
[22| 130]). but they appear not to have been widely used for models in archaeology. 
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See Fig. [2]for the Aegean when D = 100km and 110km. [The reader should impose 
the pattern of sites in Fig. [2] on the geographic template of Fig. \T\ to identify 
the sites.] At this value the S. Aegean splits into four identifiable zones; Crete, 
the Cyclades, the Dodecanese and the Peloponnese. If we increase the distance 
to D = 130km, these regions connect [not shown]. These are relevant values for 
the MBA, for which the dominant marine technology is sail, with vessels routinely 
capable of travelling large distances in single journeys. As a marker we note that 
the distance from Knossos to Akrotiri is in excess of 100km, which we believe was 
achievable in a single journey. Typically we take D to be 100km or somewhat 
greater. 



o 




Figure 2: On the left an example of an MDN (Maximum distance network) while the right 
hand network is a k-nearest-neighbour graph or PPA (Proximal Point Analysis) network. 
These are for the 39 MBA Aegean sites of Fig. [Tj using estimated travel times where 
land travel has a friction coefficient of 3 compared to sea travel. In the MDN network, 
thick edges in black are present if the separation is less than D = 100km (thin grey lines 
indicate edges between 100 and 110km). In the PPA network thick black lines indicate 
the case where each vertex is connected to its three nearest neighbours k = 3 (thin grey 
lines indicate connections to fourth and fifth nearest neighbours). In both cases the colour 
of vertices indicates the connectivity of the network defined by the thick black edges only. 

The second thresholding method is more common in archaeology where it is 
known as PPA (Proximal Point Analysis) [311 HS1 US 12 132] though it has received 
less attention in other fields, e.g. as 'k-nearest-neighbour graphs' [211 H]. Here each 
site is connected to its k nearest neighbours (k is a parameter of the model) to 
give a directed network, see Fig. El though the directions are usually ignored in the 
literature. An example is given in Fig. [2] for k = 3 nearest neighbours. It differs 
strongly from its MDN counterpart in that, by definition, sites on the extrema of the 
map that are separated from their neighbours will connect nonetheless, despite the 
distances, because the assumption is that they will connect somehow, independent 
of how easy this may be. As a result, there is a tendency for sites to be connected 
in 'strings'. 
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Figure 3: A simple Gravity Model with edge weights Fy = SiSjf(dijfD) for the 39 
MBA Aegean sites of Fig. [H using estimated travel times where land travel has a friction 
coefficient of 3 compared to sea travel. The function f(x) used is given in ([1]) with a 
distance scale of D = 100km. Vertices are proportional to the site size (0.33, 0.67 or 
1.00 for small, medium and large sites). On the right all edges are shown but with their 
thickness proportional to the flow Fy . The same network is shown on the right but now all 
links with F{j < 0.23 are shown only as thin grey lines. In addition the colour of vertices 
indicates the connectivity of the network of strong (i^,- > 0.23) links. 

The MDN can be thought of as imposing prohibitive costs on implementing long 
single journeys and PPA as imposing prohibitive costs on sustaining more than a 
few social interactions. The last simple model worth mentioning is the elementary 
'Gravity model' (TQl EQ] . For this the flow Fij from site i to site j is assumed to take 
the form Fy = SiSjf(dij/D) (with Fa = 0) where the distance 'potential' f(x) is a 
monotonic decreasing function, reflecting the ease of travel from i to j or, as we shall 
see, the effective cost of travel. In the context of the MBA Aegean we want short 
trips by sea to be relatively easy, or of low cost, while we want a strong cutoff or 
penalty for single trips of distance D or more. We stress that, unlike for the MDN, 
D is no longer an absolute cut-off, but, with /(l) = 0.5, sets the scale above which 
single journeys become increasingly difficult. D therefore takes somewhat smaller 
values than we would read off from the MDN. 

We have chosen a generic form for f(x) that is a smoothed out version of the 
MDN , namely 



In our work we have used only f3\ = 4.0 and /3 2 — 1-0. As long as there is a plateau 
for x < 1 followed by rapid falloff for x > 1 the outputs are similar. As we see in 
Fig. 3, it produces a dense network, although many of the long distance links will 
be weak. In particular either of the previous two thresholding techniques could be 
used to create a network from the flows F^. In this sense this Gravity model can 
be seen as a generalisation of the previous two models to sites of different sizes. 
In the first two models sites are of fixed and equal size. This, and subsequent, 




(1) 
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Gravity Models avoid this limitation. In particular, by aggregating island commu- 
nities into an island 'centre of population' and a similar conflation of inter-island 
exchange, they assume that, at the island level, the whole is the sum of the parts. 
By this coarse-graining, they minimise our need to have detailed local knowledge, in 
that sense making the best of poor data. However, the elementary model above fails 
to capture the rich variety of possible interactions — there is no feedback between 
interactions and site size. In arriving at the flow from % to j, we could change the 
location of all other sites and it would make no difference. 

4 Optimal Geographical Networks 

The models which are most fit for our purpose are ones in which the strength of the 
interactions reflect both the local geographical topology between two sites but also 
the wider regional structure in which these sites reside. In different ways the two 
model types that have been developed most fully are optimal models, either looking 
for the 'most likely' networks, all other things being equal, or the most 'efficient' 
networks. 

The simplest example of the former is that of the doubly constrained Gravity 
Model (DCGM). Here the flow takes the same form as before but with coefficients 
determined self-consistently:- 

F l3 = AtQiBjIjffe) , 4; = E B k I k f& , -7T=E A k O k f(^) . (2) 

1 k 3 k 

Although not immediately transparent, solutions of this form are in fact turning 
points with respect to the N(N — 1) variables (as F# = 0) of the Hamiltonian 

(ij) i j j « 

-P[C-Y,(Fii*i)]. (3) 

ij 

The first term in (J3J) is the negative of the entropy of a network of flows Fy. 
Maximising this entropy term produces the most likely distribution of the total 
exchange, J^ijFij, amongst all the possible edges, in the absence of any further 
knowledge as to how exchange occurs. As all edges are treated equally by this 
entropy term, it is the remaining terms which produce a non-trivial solution. 

The iV parameters {a^} are Lagrange multipliers which enforce constraints on 
the total outflow, Oj = J2jFij, from each site i. The {ct'j} do the same job for 
the total inflow of site j, Ij = Yli^ij- Here O, and Ij are input parameters of the 
model which we have taken proportional to the carrying capacity/population. The 
simplicity of the model allows for an algebraic solution for these Lagrange multipliers 
which allows us to replace the {a, c^} with the normalisations {A} and {Bj} in (j2J). 
Enforcing these constraints ensures the flow along the edge from i to j depends 
directly on the flow from i to all other sites and upon the flow into j from all other 
sites. The largest contributions will be from sites near to i and/or to j so the 
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interactions do depend on the whole region. Changing the location of other sites in 
the region will alter the flow from i to j, the type of behaviour we are seeking. 

The last Lagrange multiplier, (3, is associated with the total 'cost' C for the 
pattern of flows. Here the cost per unit of flow is given by Cjj for a link from i 
to j. The total cost C could reflect many factors rather than monetary and so it 
is rarely known in a physical example. So rather than specify the total cost as a 
parameter of the model, it is normal to keep /3 as a parameter. Larger (smaller) (3 
means higher (lower) costs. In practice (3 dependence becomes absorbed into one 
or more parameters into the travel, or distance, 'potential' f(x) of ([T]) on making 
the identification f(dij/D) = exp(— 0Cij). Thus larger D means cheaper costs and 
corresponds to larger (3. 

An exemplary network of the DCGM is given in Fig. HI The drawback of this 
model is that we have to fix both input and output site sizes and these do not 
respond to the pattern of interactions which emerges. For example, in our MBA 
Aegean context, Knossos can never alter the volume of economic, cultural or social 
exchange from the values we gave it at the start. Further, because inflows and 
outflows are fixed, we get a repeat of the pattern in the PPA that remote sites are 
still strongly connected, even if distances are large. 




Figure 4: Two different optimal Gravity models for the 39 MBA Aegean sites of Fig. [TJ 
using estimated travel times where land travel has a friction coefficient of 3 compared to 
sea travel. In both cases the distance scale is D = 100km. Vertices are proportional to the 
site size (0.33, 0.67 or 1.00 for small, medium and large sites). On the left solutions for 
the doubly constrained Gravity model (DCGM) are shown. On the right is the Rihll and 
Wilson Gravity model (RWGM) © with 7 = 1.3. 

Some of these issues are addressed in the work of Rihll and Wilson [251 ES] who, 
instead, use a model originally devised to study the emergence of dominant retail 
centres j6j [Ml ESJ [U |9], which we now summarise. The basic Rihll and Wilson 
Gravity model (RWGM) gives the flow from site i to site j as 






(4) 
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As before, Oi is the total outflow of site i, Oi — J2j Fij an d Ij is the total inflow 
of site j, Ij = The difference from the DCGM is that, while outflows Oi 

are still input parameters of the theory (taken to be proportional to site carrying 
capacity), the inflows Ij are now outputs determined by the model. These are used 
by Rihll and Wilson to assign an importance to a sitqj. Solutions to this nonlinear 
equation are in fact turning points of the Hamiltonian 

H = y,F,,(\n( I-]/) - 1) - 5>[o, - - 0[C - Y,(F ijCij ) 

i,j i j i,j 

- 7 [X-£i^lnOTF fcj )-l)], (5) 

i,j k 

again understood as maximising entropy. As before, we have encoded (3 and the cost 
function Cy in the distance potential / of (QQ). Likewise we retain the 7 Lagrange 
multiplier in the solution (jlj) rather than specify the unknown X in the Hamiltonian 
(J5]). An important solution is where the flow into one or more sites is zero so that 
Ij = Fij = for a finite number of sites j, but for all values of index i. 

In principle there are many ways to find these solutions. Rihll and Wilson 
generate a sequence of values {lj(0), Ij(l), . . . , Ij(t), . . .} where 

Ij(t + 1) = Y.Mt)O l (I J (t))y(d lJ /D) , A-\t) =Y.Il(t)f{d lk /D) (6) 

i k 

Provided the parameters are chosen suitably the limiting value of this sequence, 
linit-j.oo Ij(t), generates flows which optimise H of Rihll and Wilson use an 
"egalitarian hypothesis" and set the fixed outputs and the initial inputs all equal 
Oi — Ij(t — 0) = 1. This is deemed equivalent to the assumption that "all sites were 
approximately equal in size and importance at the beginning of the period under 
consideration", since there is often little information available to make any other 
hypothesis. 

In principle, t could just represent some computational time and have no physical 
significance. However it is sometimes interpreted as physical time |34j in which case 
this represents a whole further set of assumptions. In particular one can view the 
finite difference equation as the Euler method for solving the differential equation 

^ = e(j2A i0i (lM J e^-KI^ (7) 

where K and e are new constants. 

An example of the RWGM is given in Fig. H]for 7 > 1, which generates hubs. 
Taking values 7 < 1 recreates networks closer in structure to those of the DCGM 
(not shown). As with the DCGM, the enforced outflows ensure that however remote 
a site may be, it will always be connected. As a result, although we have not been 
able to show this in the limited Figures in the text, these models give networks 
whose important links are largely insensitive to our choice of D, within reason (e.g. 
varying D from 80km to 130km). 

Whatever choice we make, one possibility is to exclude such peripheral sites, as 
we have for our chosen MBA sites and for the Iron Age Greek mainland sites chosen 
by Rihll and Wilson. 



1 In a second self-consistent version, Rihll and Wilson set Oi — Ii- The model parameters now 
include an initial value for for each site. We will not consider this variant further. 
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5 ariadne 



We conclude with a discussion of the model, named ariadne, that we have developed 
over a period of time for examining MBA maritime networks in the S. Aegean 
[121 [T7]. We wanted site sizes to be informed, but not absolutely constrained, by the 
local resources, with physical geography placing similar limitations on interactions. 
Rather than look for the most likely networks, as we have been doing with gravity 
models, we are looking for the most efficient networks from the viewpoint of the 
costs and benefits that the network demands and provides. 

In our model the input parameters are the fixed local site resources, Si, and 
the site locations plus four independent parameters characterising the component 
parts of the 'Hamiltonian' or 'social potential' H, understood as representing the 
costs minus the benefits in sustaining the network. As outputs the model produces 
the flows Fij between different sites i and j but now also produces self-loop terms 
Fa > to represent parts of the population not involved in interactions. The values 
are found by minimising the following HamiltoniarJl H (i.e. maximising benefits 
with respect to costs):- 

H = -KY^m(l-W i /S i )-Xj^Fi j f(d ij /D)W :i + Jj2Wi + fi^2F ij . (8) 

i i,j i i,j 

where Wi = J2j Fij represents a variable total population at site i. Note that in this 
model Wi is not the same as the total output Oj = J2j\j^i Fij — Wi — Fa if some of 
the population at site i stays at home, Fa > 0. 

The first term of H is proportional to k is a measure of the benefit of local 
productivity. It is minimised by setting Wi = Si/ 2, only over-exploitation Wi > 
Si produces positive energy (negative benefits). This term is similar in form to 
Wj[ln(Wi) — 1] which would be present if we were maximising the entropy associated 
with the arrangement of the total 'population', J2j Wj, amongst the sites. This is in 
contrast to the optimal Gravity models where it is the entropy associated with the 
distribution of people amongst the edges which is maximised. 

The second term, proportional to A, provides benefits of exchange, as tempered 
by geography through the function / of (JTJ. Links over distances much longer than 
the parameter D produce relatively little benefit to reflect the low probability of 
such direct links being maintained. The benefit of a link is also deemed to be in 
proportion to the target site's total size Wj. As the source site's size is related to 
the flow from i to j, F^, this term is reminiscent of the product of source and target 
sites sizes that characterise gravity models as in (J2J). However our product is in the 
Hamiltonian, not in the solution for the flows. Our nonlinear dependence on the site 
sizes in this term is a key distinction between our models and the optimised gravity 
models presented earlier. Thus, an increase in interactions from i to j will produce 
a positive feedback and increasing site j's size will produce even greater benefits, 
encouraging in turn more growth in the flow between the two sites, possibly (but not 
necessarily) requiring an increase in the size of site i. The first term proportional to k 

2 In previous papers [HJ[I7J[S7] we worked in terms of fractional values for total site size Vi and 
for flows ey with Wi = SiVi and Fij — Sii^e^-. The Hamiltonian is then H = —k'Y)- 4SjVi(l — 

v i) ~ ^J2i,j( S i V i) e ij f( d ij / D )( S 3 V j) + J J2i S i V i + M SiVitij. 
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prevents this process from running away but it does allow the volume of interactions 
to grow if the dynamics favours itQ 

The last two terms, proportional to J and fi are the costs in sustaining the 
population and the total flow. The J term is equivalent to setting the total output 
rather than individual site outputs in Gravity models, so that J = i n the 

optimal Gravity models of ([3]) and (J5]). The term proportional to \x is different only 
if Fa 7^ for some sites i, emphasising that the full potential for interactions need 
not be met in our model — some 'people' can stay at home and need not 'trade' 
if there is insufficient benefit. This is another improvement over gravity models as 
in ariadne a remote site i separated by several D from the nearest neighbour will 
probably not interact, Fij = if i ^ j, and will probably maintain a population of 
about Si/2 based on the benefits of local resources encoded through the k term. 

We also impose some further constraints on our variables. Clearly we demand 
that Wi > and F^ > 0. More importantly we also impose a short range cutoff 
such that, for sites separated by less than a certain minimum distance, cfy > d m i n , 
we set F^ = 0. Technically, if we impose this cutoff then we can split any site into 
two pieces separated by less than d m i n , and these two small sites behave like a single 
unified site. In this way we have an explicit scale for our coarse graining. We suggest 
that this scale be set by the distance of travel in a day by land, say 10km. In our 
case it only effects three sites on the Southern coasts of Crete. 

We note that the scale of the Hamiltonian is irrelevant so we usually choose 
k = 1.0 leaving us with the three Hamiltonian parameters A, J, \i (the last two may 
be of any sign) plus the distance scales D and d min . The final parameters are the 
fixed site resources Si along with their separations dy. However, unlike for the 
previous models, the networks are now more sensitive to the values of D chosen. 
For D noticeably less than 100km the networks becomes fragmented, whereas for D 
more than 130km, say, the networks have a tendency to become very dense. This is 
a direct reflection of the typical distance scale ~ 130A;m (not quite identified with D, 
see earlier) for which the sites form a connected whole. We can invert the problem 
to say that the existence of a strong maritime network requires a marine technology 
whose single journeys match the distances required to produce a network. As we 
noted earlier, this is plausible, given the 100km plus of the distance from Knossos 
to Akrotiri. 

To produce our networks, we find a set of values for our N 2 parameters F^ 
which produce an approximate minimum for H of (|S]) using a Monte Carlo method. 
Unlike the RWGM, ariadne is not strictly deterministic, and only discriminates 
between comparably optimal solutions statistically. This introduces a stochastic 
element which we feel is appropriate given the uncertainties in modelling such a 
complex system. Of course we do not know the values for our input parameters, so 
our approach is to look at how the networks change as we change our parameters. 
Some examples are given in Fig. [5j For instance reducing the benefits of interaction, 
reducing A, destroys the longer distance weak links that are maintaining the global 
connectivity. This also produces a general reduction in site sizes despite the main- 
tenance of more localised networks. Likewise, the variation in city sizes, a feature 

3 While one could imagine other dependencies on the target site size, say {Wj) 1 with some new 
model parameter 7, we have not found such an additional parameter necessary. 
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which emerges in this period and exemplified by the size of Knossos, can be pro- 
duced by emphasising the non-linear terms in our Hamiltonian (also an attractive 
feature of the RWGM). Finally, though we produce equilibrium networks, we can 
tackle questions about time evolution. Slow evolution can be simulated by com- 
paring different values for our parameters, mimicking adiabatic changes in physical 
thermodynamic problems. However we can also look at 'quenches', as in removing 
Akrotiri after the eruption of Thera |18j . 

6 Conclusions 

Despite the well-established nature of many of the network modelling approaches 
discussed here, a detailed statistical comparison remains to be given. The main 
purpose of this article has been to start this process by considering the relationship 
among popular archaeological models which take geography as their primary driving 
force. Some of their key properties are summarised in table [2J 



Model 




Network Type 


Site Sizes 


Type 


(Fig.) 


W D Other 


Outflows 


Inflows 


Voronoi 


none 


Fixed Equal 


XTent 


UW D Trees 


Fixed & Different 


MDN 


© 


UW UD 


Fixed Equal 


PPA (ED 


UW D fc out >0 


Fixed Equal 


Simple GM © 


W UD 


Fixed Different 


Fixed Different 


DCGM 


© 


W D k in ,k out >Q 


Fixed Different 


Fixed Different 


RWGM (UD 


W D A; out >0 


Fixed Different 


Variable 


ariadne 


© 


W D 


Variable 


Variable 



Table 2: Summary of some of the features of different models of interactions for geograph- 
ically embedded systems. (U)W = (un)weighted, (U)D = (un)directed. k > (k out > 0) 
indicates that a model forces all sites to have at least one (outgoing) edge, however isolated 
it may be. 

Our emphasis has been on the last two models, the Rihll and Wilson Gravity 
Model (RWGM) [25| 126] and our own response to these issues, the ariadne model 
[T2j [T7] . Both of these use optimisation as a key principle, drawing on the wealth of 
experience from statistical physics. However it is interesting to highlight what we 
consider to be positive aspects of our ariadne model, and how these compare with 
the other models discussed here:- 

• ariadne and the gravity models give weights to interactions so that we can 
have both strong and weak links in the nomenclature of Granovetter [H] . 

• ariadne gives the most likely arrangement of 'population' over distinct sites, 
whereas Gravity models optimise entropy of links. 

• ariadne has no absolute constraint on individual site sizes. In contrast DCGM 
constrains both inflow and outflow while RWGM does not constrain inflow at 
sites. 
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• The more isolated a site is in ariadne, the less it will participate. However in 
gravity models and PPA, the constraints force even the most remote sites to 
be integrated into the system. 

• More generally, both the RWGM and ariadne have to be understood statis- 
tically. From the viewpoint of ensemble theory RWGM is a microcanonical 
description of network flow, whereas ariadne provides a grand canonical de- 
scription. 

It can be argued that different models are capturing different types of phenom- 
ena. For instance, the XTent and RWGM (for 7 > 1) are defining zones of control, 
which sites dominate their neighbours, rather than defining interaction patterns. 
In practice, the strengths of models only become apparent in specific applications 
to the archaeological record. Whereas PPA has been used for EBA Cyclades with 
rowing technology [I] we have used ariadne to describe the MBA S.Aegean with 
its sailing technology [121 H7J |I27J IB] • The maritime networks here are characterised 
by an ability to make single journeys on a scale that is the typical length for which 
the S.Aegean is largely connected. For networks requiring travel on a scale of sev- 
eral days journey (e.g. the East Mediterranean in the LBA) the RWGM may be 
appropriate. We shall consider this elsewhere. 

There are many properties that we have not addressed. In particular, the strongly 
non-linear behaviour of RWGM and ariadne mean that the resulting networks have 
a propensity for instability. For our model these are given in more detail elsewhere 
P2J [T7J |27J [18] to which we refer the reader. 
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A Supplementary Material 

The following tables contain all the information needed to generate the examples in 
this paper. Simple replacements of text, e.g. using a text editor, should be able to 
convert the ffl^X source into other suitable formats. 
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Number 


Name 


Size 


Lat 


Long 


1 


Knossos 


1 


35.3 


25.16 


2 


Malia 


1 


35.29 


25.49 


3 


Phaistos 


1 


35.05 


24.81 


4 


Kommos 


0.67 


35.02 


24.76 


5 


A.Triadha 


1 


35.06 


24.79 


6 


P-kastro 


1 


35.2 


26.28 


7 


Zakros 


0.67 


35.1 


26.26 


8 


Gournia 


1 


35.11 


25.79 


9 


Chania 


1 


35.52 


24.02 


10 


Akrotiri 


0.67 


36.35 


25.4 


11 


Phylakopi 


0.67 


36.73 


24.42 


12 


Kastri 


0.67 


36.22 


23.06 


13 


Naxos 


1 


37.11 


25.38 


14 


Kea 


0.67 


37.67 


24.33 


15 


Karpathos 


0.33 


35.42 


27.15 


16 


Rhodes 


1 


36.42 


28.16 


17 


Kos 


0.67 


36.88 


27.28 


18 


Miletus 


1 


37.84 


27.24 


19 


Iasos 


0.67 


37.28 


27.42 


20 


Samos 


0.67 


37.66 


26.87 


21 


Petras 


1 


35.2 


26.12 


22 


Rethymno 


1 


35.35 


24.53 


23 


Paroikia 


0.67 


37.08 


25.15 


24 


Amorgos 


0.33 


36.82 


25.15 


l 25 


los 


0.33 


36.73 


25.29 


26 


Aegina 


0.67 


37.75 


23.42 


27 


Mycenae 


1 


37.73 


22.76 


28 


A.Stephanos 


1 


36.8 


22.58 


29 


Lavrion 


0.67 


37.7 


24.05 


30 


Kasos 


0.33 


35.42 


26.91 


31 


Kalymnos 


0.33 


36.98 


27.02 


32 


Myndus 


0.67 


37.05 


27.23 


33 


Cesme 


0.67 


38.32 


26.3 


34 


Akbuk 


0.67 


37.41 


27.41 


35 


Menelaion 


0.33 


37.11 


22.37 


36 


Argos 


0.67 


37.63 


22.73 


37 


Lerna 


0.67 


37.5 


22.73 


38 


Asine 


0.33 


37.56 


22.86 


39 


Eleusis 


0.67 


38.04 


23.54 



Table 3: List of the positions and sizes of the 39 MBA Aegean sites. 
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Figure 5: Exemplary networks in ariadne for the 39 MBA Aegean sites of Fig. [TJ The 
distance scale used is D = 100km with distances based on estimated travel times in which 
land travel has a friction coefficient of 3 compared to sea travel. Vertices are proportional 
to the site size (0.33, 0.67 or 1.00 for small, medium and large sites). Top right has 
J = -1.00, /J, = 0.500, k = 1.00 and A = 4.00. Top left has A = 2.50 but other values 
unchanged. This shows the effect of reducing the benefits of interaction on the network. 
The lower right network has J = -0.975, fi = 0.500, k = 0.90 and A = 4.00. This choice 
ensures top right and bottom left networks have Hamiltonians with the linear Wi coefficient 
(this is effectively —4k + J + /i as here Fa m 0) but the coefficient of the quadratic term 
in Wi is lower in the lower network encouraging larger variations in site sizes. The lower 
left network has exactly the same parameter values of the top right but site 10, Akrotiri 
on Thera, has been removed to reflect the post-eruption geography. 
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^ilG_ 


300 


4 , o 


1 i 5 


1 i 5 


2 i 9 


291 


2, [ 1 3 


2, 1 , 


2 o 


212- 


898 


152 


2 8 o 


2 ( , o 


221 


94 


37 


88 





288 


88 2 


1 3 


222. 


-12 4 - 


88- , 


2-1— 


-112. 


^4 


2-2. 


81 


76 


222- 


70 


5 2, 5 


222. 


■ '< 8 . , 


2, , 4 


3 2, 


J2i 


_U£ 


65 


222-* 


_21_5^ 


_21_5 i 


-12- 


48 


39 


_2_U 


221* 


221*: 


303 


221 




108 


221 


2212 


-2^4 


444 


2 8 : 





157 




197 


18 


37o 


— i 




88, 


76 


221 


22), 


8 ■ , 8 


2i2^ 


40 5 


221 




391 


88 8 


22 


80 


95 


135 


1 3 5 


1 3 5 


180 


195 


144 


60 


135 


1,8 


161 


2(19 


2(8 ) 


2,3 


3 5 1 


', 2 1 1 


840 


8 2 


88 2 


157 





— — 

201 


2)10 


104 


281 


341 


2 2, 5 


208 




290 


; ', ( 1 s 


808 


349 


828 


3 1 3 


— — 

291! 


280 


309 








380 


8)80 


2,80 


233 








86 


84 




13 


119 


203 


287 


197 


188 




163 









62 






310 






28!) 


188 


203 


168 


200 


375 


253 






185 


24 


-'!)!! 


183 


305 


3 5 


8, 5 




217 


o u 2 




79 


120 


264 


67 


1 75 


229 


244 


140 


148 


170 


132 


197 


, 1 } , ; 


02 





67 


2 5 4 


849 


3 1 5 


1 3 


199 


124 


141 


175 


104 


408 


''57 


1 3 ."■ 


■>■'! [ 


197 


25 


108 


100 


333 


333 


333 


195 


208 


180 


177 


47 


70 


208 


53 


145 


282 


205 


191 


215 


230 


194 


185 


104 


44 


67 





207 


302 


200 


159 


204 


177 


192 


211 


22!) 


3,5!) 


259 


249 


250 


212 


26 


325 


330 


490 


490 


490 


395 


400 


2,08 


250 


258 


148 


187 


198 


90 


433 


402 


307 


350 


380 


335 


375 


281 


210 


254 


207 





170 


227 


77 


410 


350 


371 


289 


383 


330 


104 


142 


131 


38 


27 


393 


401 


558 


558 


558 


471 


484 


447 


304 


384 


22,5 


218 


305 


228 


522, 


559 


477 


471 


501 


443 


452 


341 


310 


34!) 


302 


170 





208 


209 


493 


407 


481 


410 


492 


330 


30 


60 


75 


190 


28 


300 


312 


308 


308 


308 


395 


409 


300 


189 


272 


198 


87 


289 


247 


455 


518 


440 


459 


490 


442 


872 


285 


2.-,2 


315 


200 


237 


208 





241 


424 


429 


445 


428 


472 


93 


233 


209 


212 


240 


29 


300 


308 


405 


405 


405 


358 


364 


332 


248 


201 


112 


194 


140 


28 


397 


492 


314 


297 


582 


275 


.3.37 


20.3 


140 


292, 


159 


77 


29!) 


241 





30!) 


290 


309 


214 


314 


334 


174 


152 


141 


77 


30 


175 


130 


238 


288 


228 


61 


68 


108 


271 


109 


205 


359 


247 


347 


30 


102 


175 


247 


240 


252 


76 


220 


23!) 


199 


204 


410 


493 


424 


309 





179 


194 


.3.32 


280 


517 


408 


440 


435 


418 


31 


205 


235 


400 


400 


400 


212 


225 


240 


828 


170 


247 


370 


170 


202 


193 


133 


35 


69 


65 


81 


210 


296 


188 


124 


177 


.350 


407 


42!) 


290 


179 





19 


175 


58 


522 


420 


404 


393 


353 


32 


28:5 


251 


415 


415 


415 


220 


238 


252 


342 


184 


201 


390 


184 


281 


190 


123 


17 


59 


49 


76 


2.30 


808 


203 


141 


192 


371 


481 


445 


309 


194 


19 





175 


47 


52,8 


438 


410 


405 


366 


33 


305 


349 


547 


547 


547 


300 


375 


373 


375 


245 


228 


372 


102 


195 


3,52, 


298 


191 


140 


197 


117 


3,52, 


808 


108 


175 


211 


289 


410 


423 


214 


.3.32 


175 


175 





184 


510 


.370 


354 


343 


295 


34 


320 


293 


457 


457 


457 


20V 


281 


295 


285 


226 


275 


429 


195 


280 


251 


170 


63 


50 


41 


70 


272 


349 


200 


164 


229 


383 


492 


472 


314 


280 


58 


47 


184 





505 


459 


437 


425 


387 


35 


398 


410 


400 


400 


400 


488 


507 


458 


287 


370 


295 


185 


387 


340 


548 


011 


539 


552 


583 


535 


405 


328 


375 


408 


359 


3.30 


330 


93 


.3.34 


517 


522 


538 


516 


505 





300 


270 


827 


33!) 


36 


366 


875 


404 


404 


404 


430 


447 


413 


277 


303 


198 


184 


270 


190 


487 


528 


470 


447 


453 


407 


424 


313 


25.", 


257 


259 


104 


30 


233, 


174 


408 


420 


43,8 


.370 


459 


300 





30 


54 


220 


37 


341 


353 


.'IV! 


2,83 


383 


414 


425 


391 


254 


281 


176 


102 


248 


108 


405 


500 


448 


425 


431 


3,85 


402 


290 


231 


235 


249 


142 


60 


20!) 


152 


440 


404 


410 


354 


437 


270 


30 





15 


158 


38 


330 


341 


373 


373 


373 


404 


415 


380 


244 


270 


105 


124 


237 


157 


454 


495 


437 


414 


420 


374 


391 


280 


220 


224 


250 


131 


75 


212 


141 


4.35 


393 


405 


343 


425 


327 


54 


15 





152 


39 


347 


352 


431 


431 


431 


417 


427 


2,8-1 


282 


252 


165 


212 


202 


91 


400 


482 


305 


348 


884 


330 


.388 


.309 


185 


197 


212 


38 


190 


240 


77 


418 


353 


300 


295 


387 


339 


220 


158 


152 






Table 4: Distances between the 39 MBA Aegean sites of Table [T] used in this paper. The column and row labels refer to the site numbers in 
Table [TJ The distances are in units of km using estimated best routes for minimum travel time. Land travel has a friction coefficient of 3 
compared to sea travel. 



